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ABSTRACT 


he paper is devoted to the derivation of approximate analytical equations for the top-of-atmosphere reflection function 
of a cloudy atmosphere. These equations are based on the analytical solution of the radiative transfer equation valid for 
optically thick clouds. In particular, we consider the radiative transfer both in the gaseous absorption bands and also in 
he regions, where gaseous absorption can be neglected. The results obtained are of importance for a number of cloud 
optics problems including cloud optical and microphysical properties determination from spaceborne optical 
nstruments. 
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1. INTRODUCTION 

Modern techniques for the retrieval of cloud parameters from satellite measurements are mostly based on the look-up- 
able (LUT) approach (Nakajima and King, 1990). However, LUTs can be substituted in number of cases by the 
analytical solutions of the radiative transfer equation valid for optically thick clouds (King et al., 1997; Kokhanovsky et 
al., 2003). This simplifies and speeds up the retrieval process considerably (Kokhanovsky et al., 2003; Rozanov and 
IKokhanovsky, 2004). The aim of this work is to introduce a number of analytical equations which can be used to model 
he top-of-atmosphere satellite signal over an extended cloud field. The work is based on analytical solutions of the 
integro-differential radiative transfer equation (RTE) derived by Germogenova (1961, 1963) for the case of arbitrary 
absorbing optically thick light scattering media. We simplify solutions taking into account that light absorption by 
droplets and crystals in clouds is weak. Also we consider the radiative transfer in gaseous absorption bands. 


2. ASYMPTOTIC THEORY 
e asymptotic solutions of the radiative transfer equation valid for optically thick media were derived by 
ermogenova (1961). They can be used for arbitrary local scattering laws and levels of absorption in the medium. The 
only limitation is that the optical thickness T must be a large number. King (1987) found that the asymptotic theory is 


accurate to within 1% for T2T', where T” =1.450 with o= (1- g). Here g is the asymmetry parameter of the 


mT 
local single scattering law defined as g = =f p(8)cos Asin @d@, where p(6) is the probability of light scattering in a 
0 


given direction or the phase function. It means that + =145 for biological media having g=0.99. Fortunately, the 


alue of g is close to 0.85 for most of water clouds in visible (Kokhanovsky, 2004a). This means that zc =10 and 
asymptotic theory is applicable to most of cases of extended cloudiness. The situation is even better for crystalline 


louds. Then the value of g is close to 0.75 and T“ = 6.0. 


he asymptotic solution for the cloud reflection function can be written in the following simple form ( Germogenova, 
1961, 1963; Sobolev, 1984; Nakajima and King, 1992): 


—2k 
mle 


-7 yz K (H) K (u). (1) 
—-le 


R( LL My PT) =R, (4, Lo» 9) i 
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Here R(/,4),Q,T) is the reflection function for a plane-parallel cloud layer. The values 4 and 4 are cosines of the 


observation 2 and incidence v4, angles, respectively, and @ is the relative azimuth. The function R (4, 4o, 9,T) is 

defined as the ratio of the intensity /, of reflected light for a turbid layer illuminated in the direction specified by the 

incidence angle 7, to the value of J, for the absolutely white Lambertian screen. Clearly, we have for the absolutely 
hite Lambertian surfaces: R =1 independently of the observation angle 2 and the relative azimuth 9. 


ere is a major problem associated with Eq. (1). It requires quite complex calculations of escape functions K ( Lt), 


reflection functions for a semi-infinite layer R „ (4,4,9) and parameters k, J, m using integral equations (Minin, 


1988). Also matrix equations can be used for this purpose( Nakajima and King , 1992). 


e problem is much more simpler in visible, where the absorption of light by clouds is negligible in most of cases and 
Eq. (1) transforms to the following result (van de Hulst, 1980): 
R (4, fg, Ps) = R? (H, Mos P) -t (T) Ko (4) Ko (4o) - 2) 
Here R° (4, Lo, Q) is the reflection function of a nonabsorbing semi-infinite light scattering layer and K, (4) is the 
escape function for a nonabsorbing medium, and 
1 
Ao aoaia] 


s the diffuse transmittance of a cloud layer under the diffuse illumination. The parameter œ =~1.07 does not vary 
onsiderably for media having different microstructures (King, 1987). Also it follows within the accuracy 2% 
(Kokhanovsky, 2004b): 


K, (u) ==[1+ 24] (4) 


at u 20.2. Effectively, Eqs. (2)-(3) reduce the problem of calculation of the reflection function of a finite cloud to that 
of a semi-infinite cloud. This is of a great importance because the function R° (,4),~) depends only on the phase 


function and this dependence is rather weak (Kokhanovsky, 2004a). Therefore, corresponding LUTs in the cloud 
retrieval algorithms can be substantially reduced (King et al., 1997). Also one can use parameterizations of the function 


R° (4,4,9), which enhances the speed of retrieval even further (Kokhanovsky et al., 2003). In particular, the 
following parameterization of this function can be used(Kokhanovsky, 2004c): 


_ A+ B( i+ My) + Cu, +F (0) 
4( u+ My) 
he parameters A, B, C and the function F (0) differ for different particulate media and can be found comparing 


RE (4, Mo) (5) 


alculations using Eq. (5) with numerical solutions of the RTE (e.g., various least square minimization techniques can 
be used for this purpose). Note that @=acos(—WU,+5ss,cos@) is the scattering angle (sg=sinv,, s = sin ® ). 
Kokhanovsky (2004b) has shown that it follows for water clouds: A=3.944, B=-2.5, C=10.664, and 
(6) = p(@)—p(@) , where p(@) is the phase function. The bar means averaging with respect to the azimuth. 
Following the same procedure as described by Kokhanovsky (2004b) we have found for crystalline clouds: A=/.247, 
B=1.186, C=5.157, and F(@)= p(@). Note that for most of applications one can neglect the small contribution due 


he last term in the nominator of Eq. (5). 


It is of importance to have a simple formulation similar to that shown above for water clouds in the near —infrared (e.g., 
for wavelengths less than 2.25 44m, where many modern spectrometers and radiometers onboard various satellite 


platforms operate). This could be done either using parameterizations of functions and parameters given in Eq. (1) 
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against the single scattering albedo @, (King and Harshvardhan, 1986) or applying so-called exponential approximation 


Zege et al., 1991). The next section is devoted to the derivation of the analytical expression for the reflection function 
of cloud fields in the modified exponential approximation. 


3. THE EXPONENTIAL APPROXIMATION 
e idea of the exponential approximation is quite simple. The main parameters and functions in Eq. (1) have following 
analytical forms for a small probability of photon absorption (PPA) #=1-—@,_ in the single photon-particle 
3(1-g) 


| a 5) 
\3(1-g) 
_ _ | p 7 
k=,3(1-g)8, l=1-44 A m=8 an (8) 


ese approximations are valid only for values of J = 0 p, /O,„ < 0.0001. Here o, and ø, are absorption and 


ext — 


nteraction event (van de Hulst, 1980): 


Fx (u)K(u), ©) 


-4 
1-2 


R, (4, llo; P) = R? (4, Lo: 9) 
) 


K(u)=K,(4 g 


extinction coefficients, respectively. For water clouds in the near-infrared the PPA # can reach 0.7 and even larger 


alues (Kokhanovsky, 2004a). So we need to consider next terms in the expansions (6)-(8) (Minin, 1991). However, 
orresponding expressions appear to be extremely complicated for most of practical applications. Therefore, Zege et al. 


1991) proposed following exponential forms for the reflection function R, (4,4,9) and the combination 
f (ilo) = mK (u)K (4) in Eq. (1): 
R., (4, Hy) = R? exp (-yu (4, Ho.) » 
f (4, 4o )=(1-exp(-2y)) Ky (4) Ko (Mo) » 


B 
=4 
N 


Ro (Mos) 
Similar exponential expression can be used for the parameter /. Namely, we have: 
l=exp(- æy ). 


u ( 4, Ly,P) = 


However, they allow us to extend the applicability of Eqs. (6)-(8) to larger values of Ø. Note that empirical 
exponential forms appear here not by chance but have deep roots related to the light diffusion theory (Rozenberg, 1961, 
Kokhanovsky, 2004b). The substitution of Eqs. (9), (10), (13) into Eq. (1) gives: 
R (H, Ho, 9, T) = Ro exp(—yu (H, Mo.) te Ky (4) Ko (4o), 
here we introduced a new parameter x= kT. The global transmittance f is given by: 
_ sinh y 
7 sinh (æy+ x) ` 
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Fig.l Dependence of the reflection on the optical thickness for wavelengths 865nm and 2130nm at the nadir 
observation and the solar zenith angle 60 degrees. The effective radius of droplets is 6 um and the coefficient of 


ariance of the gamma particle size distribution is 1/7 . Lines give the results according to the MEA. Points show 
results of exact calculations using the vector radiative transfer code SCIAPOL. The following values of the refractive 
index m of water droplets have been used: m=1.324( A = 864nm) and m=1.29-0.0004i( A = 2130nm). 


Eq. (14) transforms into Eq. (2) (and also Eq. (15) transforms into Eq. (3)) as #=0. However, Eq. (14) unlike Eq. 
2) allows to consider absorbing media as well. It is important that no new angular functions arise in Eq. (14) as 
ompared to Eq. (2). This is in contrast with Eq. (1), where parameters and functions have an implicit and complex 
dependence on the PPA £ . Eq. (14) can be used for the rapid estimations of light reflection from cloudy media and 


also for speeding up cloud retrieval algorithms (Kokhanovsky et al., 2003). The range of applicability of the 
exponential approximation (14) can be extended using correction terms derived from the numerical solution of the 
radiative transfer equation. In particular, we find that the accuracy of Eq. (14) for cloudy media can be increased using 
following substitutions: u > u(l—0.05y) , t > t—A, where 


q — Sb Melo + CH Ho 
T? 
and a=4.86, b=-13.08, c=12.76. Therefore, the final equation can be written as 
R (H, My, 9, T) = R? exp(-y0—0.05y)u ( 4, ly.) - t-A) ™ K, (4) Ko (Mo) - 
Eq. (17) is called the modified exponential approximation (MEA). We show the accuracy of the MEA given by Eq. 
17) with account for Eqs. (4), (12), (16) in Figs.1, 2 for the nadir observation conditions, the solar zenith angle 60° 
and wavelengths 865nm and 2130nm. These wavelengths are often used in cloud retrieval techniques. Note that the 
single scattering albedo is equal to 1.0 and 0.9872 at these wavelengths, respectively. The asymmetry parameter is 
0.8435 for the smaller wavelength. It is 0.8054 for the wavelength 2130nm. Exact data shown in Fig.1 are obtained 
sing the vector radiative transfer code SCIAPOL developed by us. The SCIAPOL code is based on the discrete 
ordinate approach and thoroughly tested against tabular results presented by Siewert(2000). It follows that the 
accuracy of the approximation is better than 6% for the cloud optical thickness T24 in the case considered. 
alculations performed for other angles show that the accuracy only weakly depends on the geometry providing that 
grazing observation and illumination conditions are excluded (Kokhanovsky and Rozanov, 2003). It means that the 
op-of-atmosphere reflectance over cloudy scenes can be accurately modeled in the framework of the MEA (even as 


exp(x) 
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optical thickness 


Fig.2. The errors of the MEA found from data given in Fig. 1. The solid line with the higher value of the single 
scattering albedo corresponds to the wavelength 865nm. The broken line corresponds to results obtained for the 
avelength 2130nm. 


ompared to the vector radiative transfer model). One can see from Fig.2 that the accuracy of the MEA could be 
ncreased if the exact result for the reflection function of a semi-infinite layer is used in calculations. Note that we used 
in Eq. (14) the following simple formula instead of Eq. (5) valid for the nadir observation conditions only 
Kokhanovsky, 2002) : 


0.37 +1.94 
R? (ui p) = 2A. , 
1+ fy 


e results obtained with this equation are close to those derived from more complex Eq. (5) at “=1, although Eq. 


(18) 


18) is a lot simpler. Its accuracy can be further increased by adding the function F =0.25p(1—arccos(, )) to the 
numerator. 


3. THE RADIATIVE TRANSFER IN GASEOUS ABSORPTION BANDS 


e exponential approximation presented above can be easily extended to account for the gaseous absorption. Then one 


should use the following substitutions in equations given above: T>T+T,, B (Oy, +O as. )/ (o + Ops ¢ ) , where 


ext 
he subscript “g” relates the correspondent value to the gaseous absorption process. The phase function does not need 
o be modified because we ignore molecular scattering. This could be easily accounted for if necessary. However, we 


account for the additional light absorption in the atmosphere above a cloud. Therefore, it follows for the cloud 
reflection function R in the gaseous absorption band : R =T, RT, , where we omitted arguments for the simplicity. The 


alue of R is given by Eq. (17) and T, =exp(—m,7,,,), j=1,2, where m =1/ 44, m, =1/4, and 


N 32 
Tans = > | Cansa (z)¢; (z)dz kd (19) 
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Fig.3 Dependence of the cloud reflection function on the wavelength in the oxygen A-band for cloud top heights equal 
o 0.5, 3.5, 5.5, 7.5, 9.5, 11.5 km at the cloud optical thickness 7T =10. The cloud geometrical thickness is equal to 

250m. The droplet size distribution and the illumination/viewing conditions coincide with that used in calculations 

presented in Fig.1. The atmospheric model used is identical to that described by Kokhanovsky and Rozanov(2004). 
ymbols give exact results obtained with SCIATRAN (Rozanov et al., 2005). Lines are plotted using Eq. (20). 


here C 


abs,i 


is the ith gas absorption cross section, N is the total number of gases present and c (z) is the 


oncentration of the ith gas at a given height. The integration extends from the upper cloud boundary position z, to the 
height of the optical instrument z,. The accuracy of the MEA for the gaseous absorption band can be increased if the 


single scattering contribution in the signal from the atmospheric layer above the cloud R, is also taken into account. 
hen it follows: 
R=TRT, +R. (20) 


he expression for R, is presented elsewhere (Kokhanovsky and Rozanov, 2004). We checked the accuracy of Eq. 


20) with account for Eqs. (4), (5), (12), (16), (17) by performing exact calculations using the radiative transfer code 
CIATRAN (Rozanov et al., 2005) for the oxygen absorption A-band located at the wavelengths 758-768nm. The 
atmospheric model used in calculations coincides with that described by Kokhanovsky and Rozanov (2004). The values 
R are averaged with respect to the Gaussian instrument response function with the half-width 0.225nm. The 
absorption by the oxygen was accounted for by using the HITRAN 2000 (Rothman et al., 2003) database in 
onjunction with the correlated k-distribution approximation (Kokhanovsky and Rozanov, 2004). To increase the 
accuracy of the model, we accounted for light scattering and absorption below the cloud layer using the approximate 
echnique developed by Kokhanovsky and Rozanov (2004). Results are given in Figs. 3-6. 


accounts for the cloud — upper atmospheric layer interaction in a first coarse approximation only (Kokhanovsky and 
ozanov, 2004). This interaction becomes more important for lower thick clouds (see Fig.6). 
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The same as in Fig.3 except at T 
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Fig. 6. The errors of Eq. (20) calculated using data shown in Fig.5 for various cloud top height positions and T = 50. 


4. CONCLUSION 

e present here the modified exponential approximation for the reflection function of a cloudy field. This 
approximation has the accuracy better than 5-10% both inside and outside gaseous absorption bands for cloudy media 
ith the optical thickness 7 >5 in the visible and near-infrared spectral regions for most of viewing and illumination 
onditions. Only the case of black underlying surface is considered. However, results can be easily extended for the 
ase of arbitrary underlying Lambertian surfaces (Kokhanovsky and Rozanov, 2004). The exponential-type 
approximations can be developed for other cloud characteristics including cloud transmission functions and 
polarization characteristics (Kokhanovsky, 2003). Similar approach could be used for studies of light fields inside 
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louds. The SCITRAN 2.0 software package used in this work (Rozanov et al., 2005) is freely available for non- 
ommercial use at the website www.iup.physik.uni-bremen.de/sciatran. This package also includes the possibility to 
alculate the cloud reflection function in the modified exponential approximation presented here. 
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